Bright High z Snla: A Challenge for LCDM? 



L. Perivolaropoulos a and A. Shaneloo f,,c 
a Department of Physics, University of Ioannina, Greece 
b Department of Physics, University of Oxford, 1 Keble Road, Oxford, 0X1 3NP, UK 
c BIPAC, University of Oxford, Denys Wilkinson Building, 1 Keble Road, Oxford 0X1 3RH, UK 

(Dated: June 1, 2009) 

It has recently been pointed out by Kowalski et. al. (arXiv:0804.4142 ) that there is 'an unexpected 
brightness of the Snla data at z > V. We quantify this statement by constructing a new statistic 
which is applicable directly on the Type la Supernova (Snla) distance moduli. This statistic is 
designed to pick up systematic brightness trends of Snla datapoints with respect to a best fit 
cosmological model at high redshifts. It is based on binning the normalized differences between the 
Snla distance moduli and the corresponding best fit values in the context of a specific cosmological 
model (eg AC DM). These differences are normalized by the standard errors of the observed distance 
moduli. We then focus on the highest redshift bin and extend its size towards lower redshifts until 
the Binned Normalized Difference (BND) changes sign (crosses 0) at a redshift z c (bin size N c ). 
The bin size N c of this crossing (the statistical variable) is then compared with the corresponding 
crossing bin size N mc for Monte Carlo data realizations based on the best fit model. We find that 
the crossing bin size jV c obtained from the Union08 and Gold06 data with respect to the best fit 
ACDM model is anomalously large compared to N mc of the corresponding Monte Carlo datasets 
obtained from the best fit ACDM in each case. In particular, only 2.2% of the Monte Carlo ACDM 
datasets are consistent with the Gold06 value of N c while the corresponding probability for the 
Union08 value of N c is 5.3%. Thus, according to this statistic, the probability that the high redshift 
brightness bias of the Union08 and Gold06 datasets is realized in the context of a (wo, Wi) = ( — 1,0) 
model (ACDM cosmology) is less than 6%. The corresponding realization probability in the context 
of a (uio,w\) = (—1.4,2) model is more than 30% for both the Union08 and the Gold06 datasets 
indicating a much better consistency for this model with respect to the BND statistic. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 

The discovery of the accelerating expansion of the uni- 
verse about a decade ago [l| has led to an intensive pur- 
sue of the physical origin of this acceleration. This pursue 
has been taking place in both the observational and the 
theoretical aspects of the problem. 

On the theoretical aspect, there has been significant 
progress made by pointing out several models that may 
produce the observed accelerating expansion and clarify- 
ing the limits of their predictions with respect to the ob- 
served expansion rate as a function of redshift. For exam- 
ple it has been pointed out that theoretical models based 
on modifications of general relativity 0] , interacting dark 
energy '<| or higher dimensional brane world models Q 
can easily predict an effective dark energy equation of 
state w(z) that crosses the Phantom Divide Line (PDL) 
w = —1. On the other hand, models based on general 
relativity that are free from instabilities [5| and conserve 
energy and momentum of dark energy have a w(z) con- 
fined in the range w(z) > —1. 

On the observational aspect there has been significant 
improvement of the constraints on the recent Hubble ex- 
pansion history H{z) coming from a diverse set of cos- 
mological observations. Such observations include direct 
geometrical probes (standard candles like Snla P, @, 
gamma ray bursts Q and standard rulers like the CMB 
sound horizon [1, Q) and dynamical probes (growth rate 
of cosmological perturbations [10| probed by the redshift 



distortion factor or by weak lensing [TTfl ) . 

All these observational probes are converging towards 
confirming the accelerating expansion of the universe as- 
suming the homogeneity of the universe. They have ruled 
out at several a a flat matter dominated universe (assum- 
ing a power-law form of the primordial spectrum) and 
they have produced excellent fits for the simplest cosmo- 
logical model predicting accelerating cosmic expansion. 
This model is based on the presence of the cosmological 
constant A and Cold Dark Matter (ACDM )fl2j. 

In view of the significant present and forecasted im- 
provement of relevant cosmological observations there are 
specific theoretical questions that are becoming particu- 
larly interesting. For example the question 'Is general 
relativity the correct theory on cosmological scales?' is 
particularly interesting but perhaps premature for the 
current status of observational data which still allow a 
considerable range of w(z) forms around the simplest al- 
lowed value w = — 1 corresponding to ACDM . A more 
realistic but equally important question for the current 
status of observational data is the following: 'Is ACDM 
the correct model of the accelerating expansion of the uni- 
verse?' This 'yes- no' question is more realistic because 
ACDM is a well defined model which makes clear and 
definite predictions that are easily falsifiable. On the 
other hand violations of general relativity can often be 
mimicked by (peculiar) properties of dark energy such as 
anisotropic stress or clustering [T^- 

Most approaches in answering the above question for 
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the validity of ACDM have focused on comparing ACDM 
with alternative models or parameterizations on the ba- 
sis of a bayesian analysis. Due to its simplicity and ac- 
ceptable quality of \ 2 fit, ACDM usually comes out as 
a winner in such a comparison [14| even though certain 
potential problems of the model on small [l5[ and large 
scales |l6j have recently been identified. 

An alternative approach discussed in the present study 
is to directly compare the real data with Monte Carlo 
simulations consisting of fictitious cosmological data that 
would have been obtained in the context of a ACDM 
cosmology. This comparison can be made on the basis 
of various statistics which attempt to pick up features 
of the data that can be reproduced with difficulty by a 
ACDM cosmology. 

The existence of such features is hinted by the form of 
the likelihood contours in various parameter planes con- 
taining parameter values corresponding to flat ACDM. 
For example, most Snla datasets producing likelihood 
contours in the Q\ — f2om parameter plane have the la 
contour barely intersect the line of flatness + f2om = 1 
at the lower left side of the contour [l?], EH 1 ■ Similarly, 
likelihood contours based on either Snla standard can- 
dles or standard rulers (CMB sound horizon or Baryon 
Acoustic Oscillations) and constraining the parametriza- 
tion [H 

w(z) — wo + wi- (1.1) 

_1_ "T" Z 

systematically have the point corresponding to ACDM 
(w ,wi) = (—1,0) at the lower right edge of the la 
contour while the best fit involves Wo < — 1, Wi > 
HE EH HH H2| • This feature has persisted consistently 
over the last decade and over different accelerating ex- 
pansion probes [2l[ (Snla standard candles and CMB- 
BAO standard rulers). Even though the statistical sig- 
nificance of these features when viewed individually is 
relatively low, their persistent appearance makes it likely 
that there are systematic differences between the cosmo- 
logical data and ACDM predictions. 

One such difference in the context of Snla data has 
been recently pointed out by Kowalsky et. al. [H[ where 
it was stated that there is 'an unexpected brightness of 
Snla data at z > 1'. This feature is even directly visible 
by observing the Snla distance moduli superposed with 
the best fit ACDM model (dashed line in Fig. 1) where 
most high z moduli are below the best fit ACDM curve 
(obviously the reverse happens at low redshifts to achieve 
a good fit). Notice that this bias is smaller in the con- 
text of a parametrization that crosses the PDL w = — 1 
(continuous line in Fig. 1). In the PDL crossing model 
we fix wq, wi and vary Oo m only, in order to mimic the 
ACDM number of parameters. 



1 Even though the first year SNLS data did not have this feature, 
there are preliminary indications that this feature will appear in 
the three year SNLS dataflil 
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FIG. 1: The Union08 18] distance moduli data superposed 
with the best fit ACDM model (ri 0m = 0.29) dashed line and 
with the best fit (wo,wi) — (—1.4,2) model (Qom = 0.30) 
continous line. Notice that at high redshifts z the distance 
moduli tend to be below the ACDM best fit while the trend 
is milder in the PDL crossing best fit model. 



This anomalous behavior of the data with respect to 
the ACDM best fit may be attributed to the systematic 
brightness trend of high redshift Snla with respect to the 
best fit ACDM model. It is likely that this bias of the 
Snla data with respect to ACDM best fit is also responsi- 
ble for the systematic mild preference (at la) of the Snla 
data for a w(z) crossing the w = — 1 line. The goal of the 
present paper is to study quantitatively the likelihood of 
the existence of the above described bias in the context 
of a ACDM cosmology. For this reason we use a statis- 
tic (the Binned Normalized Differences (BND)) specially 
designed to pick up systematic brightness trends of the 
Snla data with respect to a best fit cosmological model 
at high redshift. 



2. THE BINNED NORMALIZED DIFFERENCE 
(BND) STATISTIC 

The main advantages of the BND statistic discussed in 
this section may be summarized as follows: 

1. It is directly applicable on the distance moduli 
data. 

2. It is a 'yes-no' statistic for each model and it in- 
volves no comparison with alternative models or 
parametrizations. Such comparisons introduce new 
uncertainties involving the suitability of the chosen 
parametrization in the infinite dimensional func- 
tional space of possible parameterizations. For ex- 
ample a parametrization choice of w(z) = wq = 
constant leads to mo = — 1 i 0.1 [la ] favoring 
ACDM more than the CPL choice of eq. (fl~T1) 
which leads to w = -1.1 ± 0.3, w x = 1 ± 2 [lj"^ 
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FIG. 2: a: The form of the Binned Normalized Difference Q(N) obtained from the Gold06[l7| dataset assuming ACDM at 
best fit (fiom = 0.34) (red line) along with a typical form of Q(N) obtained from Monte Carlo data based on the best fit 
ACDM and the Gold06 dataset (blue dashed line), b: Similar plot testing the PDL crossing model (wq,Wi) = ( — 1.4,2) (best 
fit fiom = 0.34) instead of ACDM . The agreement between real and Monte Carlo data is significantly better in this case. 



3. It focuses on specific features of the data with re- 
spect to the best fit model (systematic brightness 
trends) thus exposing clearly the weak points of the 
model. 

4. It is insensitive to the uncertainties of the matter 
density in the sense that it does not require fix- 
ing f^om to a value motivated from other cosmo- 
logical observations or even marginalizing with re- 
spect to it to smooth out the dependence on fiom- 
In order to test LCDM (or the other similar one 
parameter parametrizations) we just minimize x 2 
with respect to fiom, find the best fit cIl(z) and 
then apply the BND statistic. No external input or 
marginalization on J7om is needed. In contrast, in 
the maximum likelihood method with multiparam- 
eter parametrizations, the choice of the matter den- 
sity, fiom can significantly affect the reconstructed 
properties of dark energy [23[ . 

The construction of our statistical analysis involves the 
following steps: 

1. Consider an H(z) parametrization (eg ACDM ) to 
be tested and a Snla distance moduli dataset /ij(zi) 
(eg Gold06[I3 or Union08[l|). 

2. Obtain the best fit form H(z) of H(z) and the cor- 
responding distance moduli best fit Jx{z). 

3. Construct the error normalized difference of the 
data from the best fit fi(z) as 



4. Consider the highest redshift bin of the normalized 
differences qi(zi) consisting of N points qi(zi) de- 
fined as the Binned Normalized Difference (BND) 

N 

Q(N)=C N J2li(zi) (2-2) 
»=i 



where Cm = —7= is a normalization factor (of no 
particular interest for our purpose), z\ is the high- 
est redshift of the dataset and the redshifts decrease 
in sequence as i increases (zn is the lowest red- 
shift of the sum). Obviously Q(N) = Q(N(z N )) 
ie the BND variable Q is also expressible in terms 
of the minimum redshift zn of the sum (|2.2p . No- 
tice that due to the central limit theorem, Q(N) is 
to a good approximation a gaussian random vari- 
able even if the luminosity distance errors are non- 
gaussian. Notice that Q(N) < implies that the 
Snla in the redshift range z\ — Zmax ^ z ^> z_/v 
are on the average brighter than the ACDM pre- 
diction. The statistical significance of this addi- 
tional brightness however requires comparison with 
Monte Carlo data (see 6 below). 

5. We increase the bin size N until Q(N) changes sign 
for the first time at N — N c , z — z c . We consider 
N c (or equivalently z c ) as our statistical variable. 
Notice that if there are systematic brightness trends 
at high redshifts then N c is expected to be anoma- 
lously large (or equivalently z c anomalously low). 

6. Finally, we ask how often does the value N c (or 
z c ) occur in Monte Carlo Snla datasets produced 
from the considered best fit cosmological H{z) 
parametrization. In particular for each Monte 
Carlo dataset we find the best fit form of H(z) 
and follow the above steps in order to find the cor- 
responding BND crossing redshift z mc . We then 
find the fraction of Monte Carlo datasets leading 
to z mc < z c (or equivalently N mc > N c ). This 
fraction is a representation of the probability that 
the dataset would be realized in the context of the 
particular cosmological model. The Monte Carlo 
realization of a given datapoint at redshift in the 
context of a particular cosmological model, is a ran- 
dom gaussian variable with mean value p(Zi) (the 
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FIG. 3: a: A histogram of the probability distribution of N mc obtained using Monte Carlo ACDM data (fio m = 0.34) in the 
context of the Gold06fl7| dataset. The thick green dashed line corresponds to the crossing redshift z c of the real Gold06 data, 
b: Similar histogram for the PDL crossing model (wq, Wi) = ( — 1.4, 2) (best fit fiom — 0.34) instead of ACDM . Notice that the 
crossing redshift z c corresponding to the real Gold06 data is a much more probable event in the context of this cosmological 
model. 
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FIG. 4: a: The form of the Binned Normalized Difference Q(N) obtained from the Union08 dataset assuming ACDM at best 
fit (flom = 0.29) (red line) along with a typical form of Q(N) obtained from Monte Carlo data based on the best fit ACDM 
and the Union08 dataset (blue dashed line), b: Similar plot testing the PDL crossing model (wq,wi) = ( — 1.4,2) (best fit 
Q m = 0.30) instead of ACDM . 



best fit distance modulus at the given redshift) and 
standard deviation <Ti(zi) (the standard error of the 
corresponding datapoint). 

Even though we have found that the BND statistic is 
particularly efficient in picking up systematic brightness 
trends at high redshifts it clearly does not consist a 
unique approach. It is possible to construct other more 
complicated statistics aiming at testing brightness trends 
or other features of the distance moduli data. For our 
purpose however which is to quantify the high redshift 
brightness trend of the data, the BND statistic is suffi- 
cient since it combines effectiveness with simplicity. 

3. RESULTS 

We have applied the BND crossing statistic to both 
the Gold06[T3| and the Union08[18] datasets testing two 



cosmological models: 

• ACDM (w = -1, wi = in cq. jjUjl ). 

• A phantom divide crossing model (wq = —1.4, 
wi = 2 in eq. (fTTT) ). 

For each model and each dataset we first obtain the 
best fit value of f2om and the best fit form of fi(z) as in 
[HI. We then construct (eq. (|2~Tj) ) and Q(N) (eq. 

(|2~2"1) ) and find N c , z c . Next we construct 1000 Monte 
Carlo datasets of distance moduli for each pair (model - 
dataset) and use them to obtain the corresponding values 

In Fig. 2a we show the form of Q(N) obtained from the 
Gold06[17| dataset (182 datapoints) assuming ACDM at 
best fit (f^om = 0.34) along with a typical form of Q(N) 
obtained from Monte Carlo data based on the best fit 
ACDM and the Gold06 dataset. Notice that the real 
data BND crossing redshift z c = 0.47 (N c = 105) is 
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FIG. 5: a: A histogram of the probability distribution of N mc obtained using Monte Carlo Union08 data under the assumption 
of a best fit ACDM cosmological model (fio m = 0.29). The thick green dashed line corresponds to the crossing redshift z c of the 
real Union08 data. Fig. b: Similar histogram for the PDL crossing model (wo, Wi) = ( — 1.4, 2) (best fit fio m = 0.30) instead of 
ACDM . Notice that the crossing redshift z c corresponding to the real data is a much more probable event in the context of 
this cosmological model. 



significantly lower (larger) than the corresponding red- 
shift (bin size) obtained from the particular Monte Carlo 
dataset (z mc = 1.2, N mc = 8). The fraction of Monte 
Carlo datasets that can mimic a BND crossing redshift 
similar to the real data (z mc < z c , N mc > N c ) is 2.2%. 
In Fig. 3a we show a histogram of the probability distri- 
bution of z mc obtained using Monte Carlo ACDM data. 
The thick dashed line corresponds to the crossing redshift 
z c of the real Gold06 data indicating that it is an unlikely 
event. Figs. 2b and 3b show the corresponding plots ob- 
tained with the same dataset (Gold06) but the tested 
model is the PDL crossing (wq,wi) = (—1.4,2) (best fit 
n 0m = 0.34) instead of ACDM (w ,Wx) = (-1,0). In 
this case we find z c = 1.01, N c — 15 (Fig. 2b) and the 
probability of z mc < z c (N mc > N c ) is 32% (Fig. 3b). 
Clearly, the Gold06 data are significantly more consis- 
tent with this cosmological model according to the BND 
statistic. 

A similar analysis as the one shown in Figs. 2, 3 is 
shown in Figs. 4, 5 for the case of the Union08 dataset 
(307 datapoints). The results are similar and consistent 
with those based on the Gold06 dataset. For ACDM we 
find z c = 0.53 (N c = 137) which is reproduced (z mc < z c ) 
only by 5.3% of the corresponding Monte Carlo datasets. 
In contrast, the crossing redshift z c = 1.37 (N c = 2) 
obtained with the (wo,wx) = (—1.4,2) (best fit figm = 
0.30) is consistent with all of the corresponding Monte 
Carlo datasets. 

The tension between ACDM and recent datasets in the 
context of the BND statistic has been verified by compar- 
ing high redshift data with corresponding Monte Carlo 
data generated in the context of ACDM . At lower red- 
shifts (z < 1) where predictions of different dark energy 
models for the luminosity distance tend to converge, we 
do not anticipate this tension to persist. In order to ver- 
ify this anticipation, we have repeated the analysis by 
starting the procedure of generating Q(N) from a lower 



redshift rather than the highest redshift point. We have 
chosen z s = 0.8 as the staring point and we have applied 
BND crossing statistic to the Union08 datasets testing 
the two cosmological models. For both ACDM and PDL 
models the crossing redshift of z c — 0.791 (N c = 3) is 
obtained that is consistent with all of the correspond- 
ing Monte Carlo datasets (Fig. 6). Thus, despite of the 
decrease of the data errors at lower redshifts, it is clear 
that the tension between ACDM and recent data in the 
context of the BND statistic is no longer present. 

An interesting question to address is the following: 
'How robust is the derived tension between the recent 
Snla data and ACDM if there are additional systematic 
errors in the data due to a possible Snla evolution?'. In 
order to address this question, we have added quadrat- 
ically an additional error of a sys = 0.20 to the error- 
bars of all data points in the Union08 sample (and also 
in its Monte Carlo realizations) and we have repeated 
the analysis. This amount of additional systematic error 
is comparable with the original error-bars of the data 
in Union08 sample. For ACDM we found z c — 0.76 
(N c = 66) which is reproduced (z mc < z c ) by 12.1% of 
the corresponding Monte Carlo data sets. It is clear that 
by assuming this additional systematic error, the consis- 
tency of ACDM and Union08 data is increased with re- 
spect to 5.3% consistency obtained from the actual data. 
In Fig. 7 we show the resulting histogram for this case. 
Clearly, the tension between data and ACDM is reduced 
by including the additional errors but it is still beyond 
the lcr level. 



4. CONCLUSION 

We conclude that according to the BND statistic, the 
Gold06 and Union08 datasets have probability 2.2% and 
5.3% respectively to have emerged in the context of the 
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FIG. 6: a: A histogram of the probability distribution of N mc , using z s — 0.8 as the staring point instead of using the highest 
redshift point. Result obtained using Monte Carlo Union08 data under the assumption of a best fit ACDM cosmological model 
(Jlom = 0.29). The thick green dashed line corresponds to the crossing redshift z c of the real Union08 data. Fig. b: Similar 
histogram for the PDL crossing model (wo,wi) — ( — 1.4,2) (best fit Qom = 0.30) instead of ACDM . Notice that both models 
show a proper concordance to the data at these lower redshifts. 
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FIG. 7: A histogram of the probability distribution of N mc ob- 
tained using Monte Carlo Union08 data assuming additional 
systematic error of a 3y s = 0.2 under the assumption of a 
best fit ACDM cosmological model (fi m = 0.298). The thick 
green dashed line corresponds to the crossing redshift z c of 
the real Union08 data with the assumed additional errors. 



various dark energy models for the luminosity distance 
tend to converge, we have verified that the revealed ten- 
sion is no longer present (Fig. 6). Our result indicates 
a potential challenge for ACDM cosmology and provides 
the motivation for obtaining additional Snla data at high 
redshifts z > 1. Such data, may confirm or disprove the 
anomalous high z Snla brightness which is mainly re- 
sponsible for the low probability of the high z Snla data 
in the context of ACDM. 

Clearly, the unexpected high z brightness of Snla can 
be interpreted either as a trend towards more decelera- 
tion at high z than expected in the context of ACDM 
or as a statistical fluctuation or finally as a systematic 
effect perhaps due to a mild Snla evolution at high z. 
However, in view of the fact that a similar mild trend 
for more deceleration than expected at high z is also ob- 
served in the context of standard rulers [1, [U , the 
latter two interpretations are less likely than the first. 

The data and the mathematica files used for the 
production of the figures may be downloaded from 
|http://leandros.physics. uoi.gr/bnd. zip| 



best fit ACDM cosmology but the corresponding proba- 
bilities for the PDL crossing model (wq,wi) = (—1.4,2) 
are larger than 30%. We have demonstrated that the 
identified tension between ACDM and recent data is due 
to the data-points at high redshift that seem to be sys- 
tematically brighter than the ACDM predictions. At 
lower redshifts (z < 0.8) where the predictions of the 
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